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Abstract 

The 2dF Galaxy Redshift Survey has now been completed and has mapped the three-dimensional 
distribution, and hence clustering, of galaxies in exquisite detail over an unprecedentedly 
large (~ 10^ Mpc^) volume of the local Universe. Here we highlight some of the major 
results to come from studies of clustering within the survey: galaxy correlation function and 
power spectrum analyses and the constraints they have placed on cosmological parameters; 
the luminosity functions of rich galaxy clusters, their dependence on global cluster proper- 
ties and galaxy type, and how they compare with the field; and the variation of galactic star 
formation activity with environment, both within clusters and in galaxy groups. 

1.1 Introduction 

Given the long and distinguished record the Carnegie Observatories have in the 
exploration of galaxies and determining their distances from us, it is perhaps fitting that this 
centennial celebration coincides with the emergence of the new generation of large galaxy 
redshift surveys that are now mapping the galaxy distribution over statistically representative 
(~ 10^ h'^ Mpc^) volumes of the local Universe. The 2dF Galaxy Redshift Survey (2dFGRS; 
Colless et al. 2001), carried out on the 3.9 m Anglo-Australia Telescope at Siding Spring 
Observatory, Australia, is one such survey. Unlike the Sloan Digital Sky Survey (Nichol 
2003), the 2dFGRS has already been completed, with its final observations taken on 11 
April 2002. At that point it had obtained spectra for 270,000 objects, providing redshifts for 
221,496 unique galaxies in the range 0.0 < z < 0.3, with a median of z = 0.11. 

The survey provides an almost complete sampling of the galaxy distribution down to an 
extinction-corrected limit of /jy = 19.45 mag over 1800 square degrees of sky. This sky 
coverage is contained within two declination strips covering 75° x 10° and 80° x 15° in the 
NGP and SGP regions, respectively, plus 99 "random" 2-degree fields in the SGR The con- 
tiguous coverage of the sky within the two strips, together with the almost complete (~ 95%) 
sampling of the galaxy distribution, has produced the highest fidelity three-dimensional (3D) 
maps ever seen of the galaxy distribution and its large-scale structure. This is clearly re- 
vealed in the cone diagram shown in Figure 1.1, where the full range of structures (knots, 
filaments, voids) are resolved in fine and delicate detail. 

One of the most conspicuous features of the galaxy distribution seen in Figure 1.1 is 
the clustering of galaxies, and the abundance of dense, "knotty" structures in which rich 
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Fig. 1.1. A cone diagram based on a 3-degree slice taken through the NGP (left) and SGP 
(right) strips of the 2dFGRS. 

clusters are embedded. It is the quantification of this clustering and the rich wealth of data 
on clusters (and their environs) that the 2dFGRS has provided that is the focus of this paper. 
Firstly, we briefly review the key findings to come from the precision measurements of 
galaxy clustering afforded by 2dFGRS, many of which have become the "flagship" results 
of the survey. We then describe the work that we have been doing on known rich clusters 
within the survey, using them as sites to study the galaxy luminosity function (LF) and 
galactic star formation in the densest environments, and hence via contrast with these same 
properties measured in lower-density regions within the survey, to draw conclusions as to 
their environmental dependence. Finally, we briefly mention work in progress on using 
the survey itself to generate a new 3D-selected catalog of galaxy groups and clusters using 
automated and objective group-finding algorithms. 

1.2 Galaxy Clustering: Key Results 

In quantifying and characterizing the clustering seen in the 3D galaxy distribution, 
2dFGRS has realized significant advances, not just in the precision of the measurements, 
but also in extending them to much larger scales. The key to the latter is the much larger 
volumes that are (sparsely) probed by the random fields (cf . the strips on their own), allowing 
structure to be measured on scales up to 400 Mpc. 

The clustering over these scales has been measured statistically using both two-point cor- 
relation function and power spectrum analyses. Being based on redshift information, these 
functions are, of course, derived in terms of the redshift-space rather than the real-space posi- 
tions of galaxies. Nonetheless, when used jointly and in combination with the matter power 
spectrum provided by cosmic microwave background measurements, they have yielded a 
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number of new fundamental results on the origin of large-scale structure, the matter content 
and density of the Universe, and galaxy biasing; these can be summarized as follows (for 
further details, see the specific references quoted and also CoUess 2003): 

(1) Unambiguous detection of colierent collapse on large scales — manifested by a flattening 
of the 2D correlation function in the line-of-sight direction at scales of 20-40 ft"' Mpc — 
confirming structures grow via gravitational instability (Peacock et al. 2001). 

(2) The detection in the power spectrum of "acoustic" oscillations due to baryon-photon cou- 
pling in the early Universe. Derivation from these oscillations of the baryon fraction: f2b/f2m = 
0. 17 ± 0.06 (Percival et al. 2001). 

(3) Measurement of 0.^ from both the power spectrum and the redshift-space distortions seen in 
the 2D correlation function: f2m = 0.30 ± 0.06. 

(4) Through comparison of the 2dFGRS power spectrum and the cosmic microwave background 
power spectrum, the first measurement of the galaxy bias parameter b* = 0.96 ± 0.08 (Lahav 
et al. 2002; see also Verde et al. 2002) and its variation with galaxy luminosity (b/b* = 
0.85+0.15L/L*; Norberg et al. 2001) and type (Madgwick et al. 2003). 

(5) Placement of a stronger limit on the neutrino fraction, < 0.13, implying a limit on 
the mass of all neutrino species of niu < 1.8 eV (Elgar0y et al. 2002). 

1.3 Cluster Luminosity Functions 

The 2dFGRS includes an abundance of rich clusters, with 947 known clusters from 
the Abell (Abell 1958; Abell, Corwin, & Olowin 1989), APM (Dalton et al. 1997) and 
EDCC (Lumsden et al. 1992) catalogs identified and further characterized by De Propris 
et al. (2002) at a time when the survey was just half complete. Of these, 60 have since 
been used for a detailed study of the cluster LF (De Propris et al. 2003). Here the selection 
was restricted to clusters with z < 0.11 in order to sample well below the predicted M*. 
A further criterion was that clusters must contain at least 40 confirmed members within 
the Abell radius (1.5/z ' Mpc). Cluster membership had been determined previously in the 
De Propris et al. (2002) study using a "gapping" algorithm to isolate cluster galaxies in 
redshift space. The clusters were also chosen to ensure a range in velocity dispersion (and 
hence mass), richness, Bautz-Morgan (B-M) type, and structural morphology. 

This sample was used to produce a series of "composite" LPs, with clusters and their 
member galaxies being combined to allow meaningful and statistically robust comparisons 
based on both global cluster properties and galaxy type. As a starting point, a composite 
LF was derived for the entire cluster sample containing 4186 members; this is shown in 
Figure 1.2. As can be seen, it provides a very high-quality "overall" cluster LF covering 
7 .5 magnitudes in luminosity (-22 < Mhj < - 1 6) . It is well fit by a Schechter (1976) function 
with a characteristic magnitude of M^^ = -20.07 ± 0.07 and a faint-end slope of a = -1 .28 ± 
0.03 (right-hand panel of Fig. 1.2). 

A series of composite LFs based on global cluster properties were also constructed, 
with the cluster sample divided into subsamples according to high (a^ > 800kms"')/low 
(ay < 800kms"i) velocity dispersion, rich/poor, "early" B-M (Types I, I-II, n)/"late" B-M 
(Types 11-111, 111), and with substructure/without substructure (see De Propris et al. 2003 
for further details). Surprisingly, no statistically significant LF variation was seen between 
these different subsamples! The only conspicuous difference seen was between composite 
LFs formed for the inner (core) and outer regions of clusters, with the former having many 
more very bright galaxies than the latter. 
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Fig. 1 .2. Left panel: The composite LF measured for 60 known rich clusters in the 2dFGRS. 
Right panel: Error contours for the M* and a parameter values from a Schechter function 
fit to this composite LF. 



Comparison of these composite cluster LFs, in particular the "overall" function shown 
in Figure 1.2, with their 2dFGRS field counterparts is readily available from the work of 
Madgwick et al. (2002). A Schechter function fit to their composite field LF over the 
same absolute magnitude range yields parameter values of M^^ = -19.79 ±0.07 mag and 
a = -1 . 19 ± 0.03. Hence, taken at face value, it would appear that the cluster LF is 0.3 mag 
brighter in M^^ and has a steeper faint-end slope (adus-ctfid = -0.1). 

However, clusters contain a very different morphological mix of galaxies compared to 
the low-density field (Oemler 1974; Dressier 1980), and before this difference in the LFs 
between the two can be interpreted as evidence for an environmental dependence, underly- 
ing LF-galaxy type effects (which have already been seen in the field; Folkes et al. 1999; 
Madgwick et al. 2002) need to be first investigated. This has been done on the basis of galaxy 
spectral type, paralleling the approach taken by Madgwick et al. (2002) for the field. Here, a 
galaxy's spectrum is typed on the basis of the relative strength of its first two principal com- 
ponents (from a principal component analysis analysis; Folkes et al. 1999), which represent 
the emission and absorption components within the spectrum. This is parameterized in terms 
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Fig. 1.3. Upper panel: Distribution of different spectral types (as measured by the 77 pa- 
rameter; see text for details) within our clusters {solid line) and in the general field {dotted 
line). Lower panel: Relationship between galaxy morphology, represented by T type, and 77. 



of the quantity 77, which is the linear combination of these two components: 77 = a pc\-pc2- 
As might be expected, clusters and the field have quite a different mix of galaxies in terms 
of 77, as can be seen in the upper panel of Figure 1.3. In line with the morphology-density 
relation (Dressier 1980) and the relationship between 7/ and morphology (T type) that is 
seen in the bottom panel of Figure 1.3, clusters are dominated by galaxies with the lowest 
?7 values (absorption-line dominated, no emission), whereas the field contains a much larger 
proportion of galaxies with higher (77 > 0) values, indicative of line emission in what are 
most likely late-type spirals. For the purposes of deriving LFs for different spectral types, 
Madgwick et al. (2002) divided the 77 scale into 4 intervals, with Type 1 galaxies being those 
in the range -5 < 77 < -1.3, Type 2 galaxies -1.3 < 77 < 1.1, Type 3 galaxies 1.1 < 77 < 3.4, 
and Type 4 galaxies 77 > 3.4. 
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Fig. 1.4. Left panel: LFs for galaxies of different spectral types within clusters. Right 
panel: The 1, 2, and 3 a error contours for the M^^ and a values from the Schechter function 
fits to these LFs. 

Using these same definitions, De Propris et al. (2003) derived composite LFs for the 
Type 1, Type 2, and Types 3+4 galaxies within their clusters; these are plotted in the left- 
hand panel of Figure 1.4. Here we see the same general trends that are seen in the field 
(Madgwick et al. 2002) in that earlier spectral types have LFs with brighter characteristic 
magnitudes and flatter faint-end slopes. This can be seen quantitatively in the right-hand 
panel of Figure 1 .4. However, careful comparison of these LFs with their field counterparts 
(Fig. 1 .5) reveals some subtle, but significant, differences, at a level of detail only a large 
survey like 2dFGRS can discern. 

These differences are encapsulated in Figure 1.6, which shows the location in the a-M^^ 
plane of the Type 1, Type 2, and Type 3+4 LFs both for clusters and the field, as well as the 
"overall" LFs based on all types. Clearly the biggest difference between cluster and field is 
for the early types, with the cluster Type 1 LF having a brighter M^^ (by ~ 0.5 mag) and a 
steeper faint-end slope (a^; - a/w = -0.5) compared to the field. In contrast, the cluster and 
field LFs for the latest types (Types 3 and 4) are, to within the errors, indistinguishable. 

The type-dependent differences in the Schechter LF parameters that are seen between 
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Fig. 1.5. Comparison of cluster and field LFs. The points represent the cluster data shown 
in the previous two figures. The thick dashed Unes in each panel represent the 2dFGRS 
field data from Madgwick et al. (2002). The thin dashed lines provide representations of the 
cluster data with no completeness correction {lower) and maximum completeness correction 
{upper); see De Propris et al. (2003) for further details. 



clusters and the field in Figure 1 .6 also put the comparison of the "overall" LFs into perspec- 
tive. Although there is a significant difference between the two, it is small in comparison to 
the differences seen between the early types (Types 1 and 2). This, together with the quite 
different mix of spectral types in clusters and in the field, would indicate that the similarity 
in the overall LFs is more by accident than design, with these mix and spectral type effects 
conspiring to produce LFs in the two different environments with characteristic magnitudes 
and faint-end slopes that are quite similar. 

In order to understand these type-dependent LF differences between clusters and the field 
in the context of cluster galaxy and environment-dependent evolution, an initial attempt has 
been made to model them using a simple "closed box" approach (De Propris et al. 2003). 
This involves making the following basic assumptions: (1) cluster galaxies have evolved 
from field galaxies contained within the volume that collapses to become the cluster, (2) this 
evolution is characterized by a suppression of star formation activity and an accompanying 
change in spectral type, and (3) the number of galaxies is conserved (i.e., mergers are ne- 
glected). Under these assumptions, the type-specific field LFs are taken to be the initial LFs 
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Fig. 1.6. Error ellipses for the Schechter function fits to the cluster (thin lines) and field 
(thick lines) LFs. Dashed lines: Type 1 galaxies; dot-dashed lines: Type 2 galaxies; dotted- 
dashed lines: Types 3+4. Only the 3 a error ellipse is shown for the field data. 



within the cluster volume, with the relative normalizations set by the mix of different types 
that are observed. 

This very simple model has some success in that it reproduces the similar shape but 
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different normalization that is observed for the Type 3+4 LFs in clusters and in the field. 
Moreover, it does cause the initial field LFs of the Type 1 and 2 galaxies to evolve toward 
the forms they are observed to have in clusters. However, at a more detailed level, it fails 
in producing too many bright Type 2 cluster galaxies, too few very bright Type 1 cluster 
galaxies, and too few faint Type 1 cluster galaxies. Obvious refinements to this model (on 
which we are currently working) are to allow for luminosity-dependent fading to steepen 
the faint-end slopes of the earlier types, and to include mergers to explain the excess of very 
bright early types at the expense of their fainter counterparts. 

1.4 Star Formation Versus Environment 

The galaxy spectra obtained by the 2dFGRS are of sufficient quality not only to 
determine redshifts but also to measure spectral-line indices and derive astrophysical in- 
formation about the galaxies themselves. Moreover, the spectra extend to sufficiently red 
wavelengths (A w 8500 A) to include the redshifted Ha emission (or absorption) line, thus 
providing a reUable means of measuring the overall star formation rate (SFR) within galaxies 
(Kennicutt 1992). In this section we describe how this has been used to track star forma- 
tion activity as a function of environment, firstly from a cluster-centric point of view, and 
secondly in galaxy groups. 

1.4.1 Within Rich Clusters 

To date, the global SFR among galaxy populations has generally concentrated on 
the two extremes of galaxy environment: the low-density field and the dense cores of rich 
clusters. For the latter, attention has been very much drawn by the discovery of Butcher & 
Oemler (1978) that such systems harbored many more star-forming galaxies in the past. But 
galaxies in cluster cores comprise ony a small fraction of the stellar content of the Universe 
and may be subject to environmental effects that are peculiar to these very high-density envi- 
ronments (e.g., ram pressure stripping, galaxy "harassment," and tidal interactions; Dressier 
2003). Much more pertinent to the evolution of the general galaxy population is the environ- 
ment between cluster cores and the field, spanning 3 orders of magnitude in galaxy density 
and which remains largely unchartered in terms of tracking star formation. 

As a first step toward redressing this situation, we have analyzed 17 z » 0.05-0.1 clus- 
ters from the De Propris et al. (2002) study, using them to trace the global galactic SFR 
continuously from their centers out to arbitrarily large radii. These clusters were chosen 
so that roughly half (10) had "high" velocity dispersions {a^, > 800kms~') and the other 
half had "low" velocity dispersions (400 < < 800km s"'). All the 2dFGRS galaxies 
within a projected distance of 20 Mpc from the centers of these clusters and within the range 
0.06 < z < 0. 10 were selected for analysis. The restriction in redshift was to limit the effects 
of aperture bias and poor sky subtraction. A luminosity cutoff of Mb, = -19 mag was then 
applied to the sample, this being the limit to which 2dFGRS is complete over the adopted 
redshift range. Finally, aU galaxies for which there were continuum/line fitting problems at 
Ha or whose EW([N ii] A6583)/EW(Ha) ratios showed evidence of a nonstellar component 
(values >0.55) were excluded from the analysis. This resulted in a final sample of 11,006 
galaxies. 

The equivalent width of the Ha line (be it in absorption or emission) was measured using 
an automated Gaussian profile fitting algorithm, which simultaneously fitted and deblended 
the Ha fine from the neighboring [N ii] A6548 A and [N ii] A6583 A lines. Since only 
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Fig. 1.7. Star formation rate as a function of cluster-centric radius, derived from our 2dF- 

GRS analysis. Filled points represent the full galaxy sample, while the triangles and crosses 
represent the "high"- and "low"-velocity dispersion clusters (see text for details). The solid 
horizontal lines represent the value of each statistic for our field galaxy sample. 



equivalent width rather than line flux measurements are possible with the 2dFGRS spectra, 
we can only infer the SFR, /i, per unit luminosity: /x/Lcont = 7.9 x 10"^^ EW(Hq:), using 
Kennicutt's (1992) SFR-EW(Ha) relation. We then normalized this to a characteristic lu- 
minosity L*: n* = fx / (Leant /i^*) = 0.087 EW(Ha), where L* is taken to be the knee in the 
luminosity function in the r' band (near rest-frame Ha), as determined by Blanton et al. 
(2001). 

In Figure 1.7, the mean and median value of /j.* is plotted as a function of the projected 
cluster-centric radius. Here the latter has been normalized by the cluster virial radius, 
which has been calculated for each individual cluster and ranges from 1.4 to 2.4 Mpc. Also 
plotted in the top panel of Figure 1.7 is the fraction of galaxies with /x* > IMq yr"\ which 
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represents the tail of the distribution, comprised of galaxies that are currently forming stars 
at a high rate relative to their luminosities. The solid horizontal line represents the values 
derived for the "field" galaxies within our sample, that is, galaxies that lie within the 20 Mpc 
selection radius, but which, from their redshifts, are identified as non-members. The brack- 
eting dashed horizontal Unes represent the 1 cr standard deviation from field to field, giving 
some estimate of the cosmological variance in the field value. 

Irrespective of which statistic is used, it is very clear that within the clusters /i* falls sig- 
nificantly below the value of the field, the difference being at a maximum at the cluster center 
and then monotonically decreasing with increasing cluster-centric radius. Importantly, con- 
vergence does not occur until R > 3/?^.' Hence, compared to the field, cluster galaxies differ 
in their mean star formation properties as far out as ~ 6 Mpc from their centers. Also of note 
is that these radial trends in /z* appear to be insensitive to cluster mass, with there being no 
discernible difference between the "high"- and "low"- velocity dispersion clusters. 

Such a radial analysis, however, is likely to be sub-optimal, since many of the clusters in 
our sample are clearly not spherically symmetric and show substructure. We have therefore 
analyzed fi* as a function of the local projected galaxy density, S, based on the distance to 
the tenth nearest neighbor (Dressier 1980). The relationship between /x* and S is shown in 
Figure 1.8, using the same three statistics that we used in the radial analysis. The vertical 
line shows the mean value of S within R^, and once again the horizontal lines show ji* and 
its 1 a variance for the field. 

We see that cluster galaxy star formation is suppressed relative to the field, the difference 
being greatest at the highest local densities, and decreasing monotonically with decreasing 
density until the two converge at S « 1.5 galaxy Mpc"^, a factor of ~ 2.5 times lower than 
the mean projected density of the cluster viriaUzed region. Yet again the trend is no different 
in the "high"- and "low"-velocity dispersion clusters, indicating that SFRs depend only on 
the local density, regardless of the large-scale structure in which they are embedded. Further 
evidence that local density is the key variant is also provided by a plot based on just those 
galaxies with R > 2Ri,, where the same trend observed for the full sample is seen. This 
would suggest that a more general view of star formation suppression be taken: that it will 
be low relative to the global average in any region where the local density exceeds a value 
of ~ 1 galMpc-2 (Mb < -19 mag). 

It is well known that galaxy morphology is very strongly correlated with local galaxy 
density (Dressier 1980), and hence consideration needs to be given to what underlying con- 
tribution this has to the ^* -S relationship seen in Figure 1.8. Using a simple model based 
on Dressler's morphology-density relation, we have calculated the expected variation in /x*; 
this is represented by the solid curves in Figure 1.8. This appears to be shallower than the 
observed relation, suggesting that the morphology-density relation is distinct from the SFR- 
density relation, with additional processes operating at S > 1 galaxy Mpc"^ that drive /x* 
significantly below that expected on the basis of the changing morphological mix. Indeed, 
it may well be the case that the suppression of star formation is the primary transforma- 
tion that occurs with envirormient, with the change in morphological mix being a secondary 
effect (Poggianti et al. 1999; Shioya et al. 2002) 

A complete description of this study and a full discussion of its implications can be found 
in Lewis et al. (2002). 
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Fig. 1.8. Star formation rate as a function of local projected galaxy density, S. The dashed 
vertical line represents the mean value of S within R^, . The curved lines indicate the variation 
expected as a result of the varying morphological mix. 



1.4.2 Within Galaxy Groups 

Extension of this analysis to the galaxy group environment has just recently be- 
come possible through the 2dFGRS being subjected to a 3D "friends-of-friends" analysis to 
identify groups (and clusters) in both position and redshift space (Eke et al. 2003). This has 
produced a catalog containing ~ 30,000 groups with at least 2 members; their distribution 
within the two 2dFGRS strips is shown in Figure 1.9. The star formation properties of this 
sample are currently being analyzed, and we show here (courtesy of Dr. M. Balogh) a cou- 
ple of very preliminary results that give the first indications of environmental trends in these 
much poorer, lower-velocity dispersion analogs to rich clusters. 

Once again the Ha line has been used as a measure of SFR, with it being identified and 
its equivalent width measured in the spectra of all group members in identical fashion to 



12 



W. J. Couch, M. M. Colless, and R. De Propris 




Fig. 1.9. A redshift slice showing the distribution of groups and clusters identified within 
the 2dFGRS using a 3D "friends-of-friends" analysis (Eke et al. 2003). The estimated 
velocity dispersion is shown by the size of the dot. 



that described above for cluster galaxies. An overall Ha equivalent width was then derived 
for each group taking the mean value of all its members. In Figure 1.10 this mean value 
is plotted as a function of group velocity dispersion, being represented by all the individual 
dots. To assess whether there is any overall trend, the data have been divided into equally 
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Fig. 1.10. The mean equivalent width of Ha for groups within the 2dFGRS versus their ve- 
locity dispersion. The dots represent individual groups whose values have been determined 
from their member galaxies; the solid line connects the mean values evaluated within equally 
populated velocity dispersion bins; the dashed horizontal line indicates the mean value for 
"field" galaxies. 



populated velocity dispersion bins and averaged; the continuous line has been drawn through 
these average values. As a comparison, the mean Ha equivalent width for "field" galaxies — 
those galaxies identified as not belonging to any group within the friends-of-friends analysis 
— is shown as the horizontal dashed line. 

It can be seen in Figure 1.10 that there is a noisy but significant trend of Ha emission (and 
hence SFR) becoming increasingly stronger with lower group velocity dispersion. Intrigu- 
ingly, however, this trend intersects the horizontal "field" line at cr « 250 km s"', indicating 
that there are not only groups whose mean EW(Ha) is below that of the field (as might be 
expected from the SFR-S relation observed in the rich cluster study), but also groups whose 
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Fig. 1.11. The mean equivalent width of Ha for binary groups with velocity dispersions 
a < 150kms~\ plotted as a function of the separation between the two galaxies. The dots 
represent individual groups; the solid line connects the mean values evaluated within equally 
populated galaxy separation bins; the dashed horizontal line indicates the mean Ha equiva- 
lent width for field galaxies. 



mean EW(Ha) is higher than the field, hinting that star formation activity is for some reason 
enhanced (cf. the field) in these lower-velocity dispersion (ct < 150km s"') groups. 

The underlying reason for this may well be contained in Figure 1.11, where we show 
the mean EW(Ha) for just the binary groups within the a < 150 km s"^ regime, plotted as 
a function of their separation (in Mpc). Yet again, the equal-numbered bin averaging and 
the benchmark "field" value are represented in the same way as in Figure 1.10. We see that 
for separations in the range 0.1-10^"^ Mpc, binary groups have, to within the fluctuations, 
the same mean EW(Ha) as the field. However, at separations below 0.1 /;"' Mpc, the mean 
EW(Ha) for binary groups is seen to be significantly higher than the field. Although this is 
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yet to be fully investigated, it is suggestive that the enhanced SFRs seen in the low-velocity 
dispersion groups (relative to the field) can be attributed to "close" binary systems where 
tidal interactions may well be responsible for this activity. 
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